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Abstract 

A relationship, linking the saturation momentum in the case of fixed and running QCD coupling, respectively, is derived from 
the Balitsky-Kovchegov equation. It relies on the linear instability of the evolution equation in the dilute regime. The leading 
orders of the saturation momenta are mapped onto each other exactly. For subleading terms a qualitative correspondence 
is achieved with a relative error going to zero for large rapidities. The relationship can also be derived for the Balitsky- 
Kovchegov equation with a cutoff accounting for low-density effects and is satisfied by the corresponding isoline functions. 
Further implications arise for the existence of travelling-wave solutions in the two situations. 



I. INTRODUCTION AND RESULTS 

In deep-inelastic scattering the total cross-section a 
for the scattering of a virtual photon with momentum q 
on a proton with momentum p is a function of the vir- 
tually Q 2 = —q 2 and the rapidity Y = ln(l/x) with 
x = Q 2 /(2p-q). At sufficiently small x < 10 -2 the cross- 
section a becomes a function of the ratio of the virtuality 
Q 2 and a function Q S 2 (Y) of rapidity Y called saturation 
momentum: a — a\Q 2 J Q 2 (Y)] pj. This phenomenon is 
called geometric scaling. Translated to the dipole scatter- 
ing amplitude N, it becomes a function of the difference 
between the logarithm L of the square of the momentum 
variable conjugate to the dipole size and the logarithm 
of the saturation momentum: N = N[L — lnQ s 2 (Y)]. 

In quantum chromodynamics (QCD) the scattering of 
dipoles is described by the Balitsky hierarchy Q. In 
the factorising limit it reduces to the Balitsky-Kovchegov 
(BK) equation [3|. In the low-density regime it in turn 
simplifies to the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
equation jj. 

In Ref. [5j scaling solutions for the BK equation are 
identified at fixed coupling and with the second or- 
der ("diffusive") approximation 01 f° r the BFKL ker- 
nel. That investigation is based on knowledge about 
the existence of travelling-wave solutions for the Fisher- 
Kolmogorov-Petrovsky-Piskunov (FKPP) equation |(|. 

This approach is not feasible for the running coupling 
BK equation because it belongs to a different universality 
class. Different methods have to be used, like a rescal- 
ing of the rapidity Y 0, introduction of curved absorp- 
tive boundaries into the BFKL equation ||, or use of a 
travelling-wave ansatz [Isj . Note that, in general, the 
" time" and " space" variables need not be linear functions 
of the rapidity Y and the momentum variable L. 

Here a mapping between the saturation momentum at 
fixed and at running coupling, respectively, is derived. 
The derivation is based on the observation that in the 
scaling regime the saturation momentum characterises 
an isoline of the amplitude N. It is derived by neglect- 
ing commutators of the BFKL kernel and the running- 
coupling function acting on the amplitude N. This pro- 



cedure is justified by the fact that the BK equation is lin- 
early unstable with respect to small perturbations around 
N = (pulled front). At this level of accuracy it is shown 
to hold for the saturation momentum in both cases at 
leading order for large rapidities Y. For the sublead- 
ing terms the qualitative features are reproduced: They 
provide a negative but small correction to the leading be- 
haviour. The relative error connected to the quantitative 
deviations vanishes for large rapidities Y. 

As said instability is already present for the BFKL 
equation the relationship also holds there. Actually the 
validity of the relation seems to be widely independent of 
the detailed form of the equation of motion and the form 
of the running coupling function's equivalent. It appears 
to hold for a larger class of differential equations. In this 
context, the linear instability seems to be a sufficient con- 
dition. The aforementioned quantitiative deviations are 
non-universal in the sense that they depend on details of 
the equation of motion which arc independent of whether 
or not it is linearly unstable. 

Furthermore, the relationship between the isolines for 
fixed and running coupling can also be derived after a cut- 
off has been introduced into the growth term of the BK 
equation in order to accomodate effects beyond the mean- 
field approximation 0, [nj, ^ . Even if, strictly speak- 
ing, the thus modified BK equation is no longer linearly 
unstable against arbitrarily small perturbations around 
N = 0, the leading-order terms for large rapidities Y 
are again mapped onto each other exactly. The sublead- 
ing terms decay exponentially |T^ . For the approxima- 
tion made in the course of the relationship's derivation 
to hold, it suffices that the BK equation with cutoff is 
unstable for sufficiently large perturbations. 

The comparison of isoline plots 1141 obtained by solv- 
ing the BK equation numerically |l5j for fixed and run- 
ning coupling respectively, according to the relations pre- 
sented below, would also provide interesting insights and 
cross checks. 

In section llll the mapping between the isolines for the 
BK equation for fixed and running coupling, respectively, 
is derived. In section llTTl the validity of the relationship is 
checked for the leading fsubsection llll A|) and sub-leading 
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terms fsubsection llll B|) of the saturation momentum. In 
subsection IIII CI the reason for the accuracy of the rela- 
tion is discussed. It is explained how it can be gener- 
alised to a larger class of differential equations. Section 
^3 treats the mapping for the BK equation with a cutoff 
taking into account low-density effects. For convenience 
appendix lAlexposes details of the connection between the 
fixed coupling BK equation and the FKPP equation. 



II. RELATION 

The BK equation for the dipole forward-scattering am- 
plitude as a function of the rapidity Y and the momen- 
tum variable L is given by: 



dN _ 
dY ~ ' 

with the BFKL kernel: 



1- 



N 



X(7)=2^(l)-^(7)-^(l~7) 



(1) 



(2) 



where ^(7) is the digamma function — the logarithmic 
derivative of the gamma function — and where the BFKL 
eigenvalue-function with a differential operator in the ar- 
gument is defined via its series expansion around 70 S 
]0;1[. 

The phenomenon of scaling manifests itself by the iso- 
lines of the amplitude N(Y, L) keeping their distance 
from each other in the direction of the variable L constant 
if the rapidity Y changes. In other words, the amplitude 
is only a function of the difference L — In Q S (Y) 2 . Conse- 
quently, in the scaling regime, the saturation momentum 
Qs = Qs(y) characterises the isoline with N[Y, Q S (Y)] = 
N s There, the equations describing two isolines dif- 
fer merely by an additive constant. 

Independent of the phenomenon of scaling, any isoline 
Li = Li(Y) of the amplitude N(Y,L) satisfies the rela- 
tion 



dU __ dN fdN 
~dY ~ ~dY [~dT, 



(3) 



For fixed QCD coupling, dN/dY can been replaced by 
the right-hand side of (JTJ. 



dY 



X 







dL 



N f - N f ' 



dN f 



dL 



(4) 



The running coupling case is obtained through the re- 
placement a — > (bL) : 



dL r 
~d~Y 



1 
bL r 



X 



d 
dL r 



N r - N r 



fdN r 



(5) 



The constant b is linked to the QCD-/3- function and 
reads: b = (lln c — 2nt)/(12n c ). n c stands for the num- 
ber of colours and nf for the number of (massless) 
flavours. 



Hence, division of Eq. by Lf, taking into account 
the common initial condition 



N f (Y ,L) = N (L) = N r (Y ,L), 



yields: 



1 dL f dL r 



LfdYo dY 



(6) 



(7) 



The constants a and b have been omitted for the sake of 
clarity. It suffices to keep in mind to exchange the two 
in a comparison. 

Analogous relations for the higher derivatives can be 
derived by repetition of the above steps: After taking 
the derivatives of Eqs. J3J and (J5J with respect to the 
rapidity Y, replace the new occurrences of the deriva- 
tive dN/dY on the right-hand side with the help of the 
BK equation Q . Subsequently, divide the expression ob- 
tained from Eq. Q by Lf. This allows to identify the re- 
spective right-hand sides at Y = Yq up to terms involving 
commutators of the BFKL kernel with the running cou- 
pling function acting on the amplitude [x(~ L _1 ]_ A. 
They are going to be omitted in what follows. Why this 
is justified is going to be investigated in subsection IIII CI 

Together with the initial condition: 

Lf(Y )=L = L r (Y ) (8) 
this hierarchy of equations can be summarised by: 
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[L f (Y )dY 



Lf(Y ) 



dY 



L r {Y ) 



(9) 



for all n G No. This defines all Taylor coefficients of 
L r (Y) at rapidity Y = Yq based on those of Lf(Y). 

Inversely, by multiplying with L r each time instead 
of dividing by Lf one obtains all Taylor coefficients of 
L r (Y) from those of L f (Y): 



dY 



Lf(Y ) 



L r (Y ) 



dY 



L r (Y ) 



(10) 



Equation can be reexpressed as: 

L r (Y + 5Y) = cxp {SY-^— -^r 
I L fK Y o) dY Q 

and Eq. I|1U|I accordingly as: 



Lf(Yo) (11) 



Lf(Y Q + SY) 



exp!sYL r (Y )4r 



( ^Y Y L r{Yv)^\L r {Y ) (12) 



for all SY . This can be verified by Taylor expansions 
around 5Y — 0. The exponentials in Eqs. i|ll|l and i|12|) 
are operators for conformal mappings. In the following 
calculations the interpretation as translation operators is 
to be used. (Alternative computations could, for exam- 
ple, involve dilatation operators.) With the definitions: 



dZf_ 
dY 



Lf(Y), 



(13) 
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and 



A. Leading order 



dZ r 



1 



dY L r {YY 
respectively, these turn into: 

L r (Y + SY) = exp^SY^L f [Y f (Z f )} = 
= L f {Y f [Z f (Y )+SY}} 

L f (Y + 6Y) = exp^5Y^-^L r [Y r (Z r )} = 
= L r {Y r [Z r (Y ) + SY}} 



(14) 



(15) 



and 



(16) 



for all SY. Yf(Zf) and Y r (Z r ) are the inverse functions 
of Zf(Y) and Z r (Y), respectively, as defined in Eqs. ljT3|) 
and (|14[) . Eqs. (|15|l and (|16fl provide a direct link be- 
tween the isolines of the solutions for the equations of 
motion in the fixed and the running coupling case. They 
hold irrespective of whether the exact or an approxima- 
tive expression, for example, the second-order expansion 
around 7 = j c , is used for the BFKL kernel 



III. COMPARISON 

For large rapidities Y, the universal expression for the 
(logarithm of the) saturation momentum — up to an ar- 
bitrary additive constant — in the case of fixed coupling 
available in the literature reads H, 0, 0] : 



Lf(Y) = 



X(7c 



-Y 



27c 



•lnF 



2tt 



7c 2 V a x"{lc) V Y 



(17) 



The BFKL kernel @ has been expanded around 70 = 7 C 
which solves the equation jcX'lrfc) = x(lc) and has the 
numerical value 7 C = 0.6275... 9]. The third term is only 
given in ^(| and there the error is of the order 0{Y r )- 
For the running coupling case the two leading terms are 
known HH: 



L r (Y) 




X"(7c) 



(is) 



with £1 = —2.338... being the rightmost zero of the Airy 
function Ai. The universally valid expression, i.e., in- 
dependent of the initial conditions, can at most include 
terms that decrease more slowly than Y^ 1 / 2 . Like in the 
case of fixed QCD coupling, this can be checked through 
a shift AF of the rapidity Y in the above equation and 
subsequent expansion around AY = 16]. Thus, in gen- 
eral, the rapidity variable Y in the last two equations 
differs by this additive constant AY. 



Starting out with the first term in Eq. (|17|) one finds 
from Eq. (|T3|) : 



ZfiX) = a 



= X(jc)y2 



27c 



Inversion yields: 



Yf(Zf) = 



27c 



(19) 



(20) 



Evaluation at Zf(Yo) + SY leads to: 



\Z,{Y„) + SY]- (2i; 



Replacing the rapidity Y in the first term of Eq. i|17|) 
by the right-hand side of the previous expression and 
afterwards a by 6 _1 yields: 

\j ( y ° + SY + Ay ) = Lr ( Yo + SY ^ ^ 

where constants have been absorbed in AY. Thus, to this 
order, Eq. (|15fl is satisfied. With analogous calculations 
also Eq. l(TB|l can be verified to leading order. 



B. Subleading terms 

With the first two terms in Eq. I|17|) one finds from 
Eq. (EU: 



Zf(Y) 



■ Xilc 

I 

27c 



■Y z 



W^ 1 



IylY)Y, 



(23) 



the exact inverse of which cannot be given analytically. 
The factor (1 — InY) varies slowly as compared to pow- 
ers of Y. Regarding it as fixed in every step, one 
can determine the inverse iteratively. Starting out with 
lnYj 0) (Z f ) = 1 leads to Eq. (231 for Yf\z s ). Replac- 
ing In Y in the logarithm in Eq. I|23|) by the preliminary 
result In Yj(Zf) [in general ]xiYj n \Zf)], subsequent in- 
version, and selection of the positive solution in every 
step finally leads to the recursion relation: 



Yj n+1) (Z f ) = - 



1 -lnYj n) (Z f ) 



2ax(7c 



A 



l-lnYj n) (Zf) 



+ [Y^(Z f )Y 



2ax(7c) 



(24) 
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It shows that the dominant behaviour for large Zt — and 

hence large Yj (Zf) — can be obtained after a finite num- 
ber of iterations. The dominant terms are given by: 



Yj n+1 \z f ) = Y} 1 \z f ) + 



2ax(7c) 



\nY^{Z } ) 



2«x(7 c ) 



2 [lnlf^/)] 2 



O 



,(25) 



where constant terms have been omitted. Replacing the 
rapidity Y in the first two terms of Eq. (|17() by the right- 
hand side of the previous expression yields: 



Lr = fi ^) (Z/) 



[\nY\ 1 \Z f )f 



7c 



8«x(7c)7c y f (1) (Z/) 



-o 



YfHZf) j 



(26) 



which can already be obtained from the second-order 
result, i.e., from Eq. (|24[l with n = 1. Evaluation at 
Zf(Yo) + 8Y leads to the replacement: 



Yj 1 \z f )^Y f ^[Z f (Y )+SY} = 



/ 2x(7c) 

«x(7c) V & 7c 



(Y + 5Y + AY), 



(27) 



which coincides with Eq. I|21|l. Finally, a has to be 
replaced by b^ 1 . 

Again, the leading term of Eq. I|18f) is reproduced. The 
first subleading terms in Eqs. l(T%|) and lf2lfl) together with 
(|27|l do not coincide exactly. However, they are similar 
qualitatively. Through their inclusion with the leading 
term, L r and L r are both diminished. The relative error 
vanishes for large rapidities Y like: (L r — L r )/(L r + L r ) ~ 
y-l/3. 

In principle, it is possible to base the above comparison 
on Eq. While the integral required for solving Eq. 

(|14f) is still known analytically for L r (Y) given by Eq. 
P&[l. the inversion of the resulting expression is more 
cumbersome than for Eq. Ij23(l . 



C. Discussion 

As demonstrated in the previous subsection, neglecting 
the commutator [%(— di), L^ 1 ]- N leads to a mapping 
exact at leading order and with subleading deviations 
whose relative error vanished for large rapidities. In what 
follows it shall be discussed why this is the case. 

Omitting said commutator is equivalent to approxi- 
mating the prefactors of the m th derivatives with respect 
to the momentum variable L occuring on the running- 
coupling side during the relationship's derivation by the 
term dominant for large L: [L^ 1 + <D{L- 2 )]d L m N w 
L~ l dL m N. Subsequently one would have to justify why 



L is effectively so large that the above steps are feasible. 
One is tempted to bring forward the fact that saturation 
physics is protected from the influence of the infrared, 
i.e., from small L. However, through the repetition of 
the identical steps the above relationship can also be de- 
rived for the BFKL equation and in the BFKL equation 
no saturation effects are encoded. 

As this line of arguments is not conclusive let us inves- 
tigate how shifting the BFKL kernel by an additive con- 
stant xil) — x(l) + ^ docs affect the expression for the 
saturation momentum. First for arbitrary values of the 
shift 8 £ M. the modified critical value 7 C for the argument 
7, defined through x(%) = %x'(%) obeys < % < 1. 
Therefore the solution stays always in the supercritical 
regime of the FKPP equation % l > 1 [2CI. i.e.. it has 
always a universal travelling- wave solution 0, |ij . Hence 
one can explore the effect of the shift 8 directly with the 
help of the expression for the saturation momentum in 
Eq. fTfl . 

For 8 > —4 In 2, xil) remains positive definite and 
the same qualitative asymptotic behaviour is retained, 
although reaching the asymptotic regime might require 
extremely large rapidities Y if x(7c) ^ 1- 

The picture changes for S = —4 In 2, where % = i and 
x(Tc) = 0, i- e ? t ne minimum of the kernel becomes zero. 
The term proportional to the rapidity Y is absent. As 
explained above, by shifting the BFKL kernel one stays 
in the supercritical regime of the FKPP equation. How- 
ever, the L and the x axes are not parallel to each other. 
In this particular situation x ~ 2t is mapped exactly onto 
L = const, (see Appendix |A"|) . whence the linear term of 
the saturation momentum vanishes although the FKPP 
equation has a supercritical travelling-wave solution. In 
other words, for 8 = —4 In 2 the BK equation is not lin- 
early unstable for small perturbations around N — even 
if the FKPP equation is. 

The effect on the expression for running coupling l|18|) 
is more drastic. While the term proportional to the 
square root of the rapidity Y vanished in unison with 
the linear term in the fixed coupling case, the prefactor 
of the second addend diverges like ~ xile) . Hence, in 
the running-coupling case the previous description breaks 
down altogether. 

Proceeding to 8 < —4 In 2, leads to ^(7) also taking 
negative values and especially to a negative critical value 
x(7c) < 0. Thus the term in the fixed-coupling satura- 
tion momentum proportional to the rapidity Y reappears 
but with a negative sign. This means that for growing 
rapidities Y the amplitude decreases. With this kernel 
the BK equation is stable against perturbations around 
N = 0. In this range, according to Eq. ^18p. the satura- 
tion momentum in the running-coupling case would even 
be complex. This shows that the latter case would have 
to be investigated anew. 

Summarising, the term of the differential equation im- 
portant for the instability around N = and hence for 
the mapping to work is the one originating from the criti- 
cal value of the kernel x(lc)- Preserving only this term in 
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the BK equation and obtaining the relationship between 
the isolines for the different couplings leads exactly to the 
previous results, because in this limit the crucial com- 
mutator vanishes [x( 7 c), L^ 1 ] N = 0. The more general 
derivation in section [H] resums additional terms which 
lead to the same overall behaviour. 

As an outlook, based on the discussions in this sub- 
section the present approach should work for any pair of 
differential equations which are linearly unstable around 
the dilute state, widely independent of the details of the 
remaining terms and the detailed structure of the equiva- 
lent of the running coupling-constant. The main step for 
adapting to another situation is replacing the right-hand 
sides of Eqs. fTljl or fT^ . respectively, by the new run- 
ning coupling function a{L) or its reciprocal [a(-L)] -1 , 
respectively. 



IV. LOW-DENSITY EFFECTS 

The BK equation describes the mean-field limit of the 
Balitsky hierarchy. The mean-field approximation is best 
satisfied in the dense regime and least in the dilute. 
There fluctuations are important. As mentioned above, 
for the relevant boundary conditions, the BK equation 
describes the propagation into a linearly unstable state. 
This leads to a high sensitivity of the solution to modifi- 
cations at the toe of the front. In Refs. it has been 
demonstrated that the principal correction to the front 
propagation speed can be simulated in a deterministic 
manner by cutting off the growth term for small values 
of the amplitude. For example in the diffusive approx- 
imation to the BK equation this leads to the following 
modified equation of motion: 



dN 
~d~Y 



X2 



d 2 N 



dN 



Xl ~dL +(^N~N 2 ) c(N) 



(28) 



with the replacement a —* (bL)^ 1 for the case of running 
coupling and where the coefficients Xii * G {1, 2, 3} are 
given in Appendix lAl and with the cutoff function: 

c(N) =9(N- a s 2 ) ■ (29) 

Note that 1 / as 2 equals the number of dipoles at satura- 
tion wherefore as 2 gives the corresponding step height. 

Carrying out the steps of the derivation beginning with 
Eq. © but this time for the modified BK equation 128|) 
instead of its standard form yields again Eqs. i|15|) 
and H16|) . Looking at the modified expressions for the 
saturation momentum one sees that Eqs. I|15|l and 
are satisfied exactly: At fixed coupling the leading term 
of the isoline equation reads \Wi fllL Il3| : 



Lf(Y) = 



Xilc) 



lcX"{lc) 



2 \n 2 {l/a 2 ) 



Y, 



(30) 



with exponentially small subleading terms [K 
ning coupling one finds |ll| : 



At 




X(lc) _ 7T 2 lcX"(lc) 

7c 2 \n 2 (l/a s 2 ) 



Y. 



(31) 



The introduction of the cutoff removes the linear insta- 
bility of the BK equation with respect to perturbations 
smaller than the threshold as 2 - However, the instability 
for perturbations larger than a s 2 is sufficient to ensure 
that the modified expressions for the saturation momen- 
tum are qualitatively similar to those of the unmodified 
version. Consequently the mapping still works. 

In this last context the term "saturation momentum" is 
avoided on purpose because after the inclusion of stochas- 
tic effects the solutions of the BK equation does no longer 
provide the observable amplitude but the amplitude for 
the scattering on a given partonic realisation of the tar- 
get [H The physical amplitude is obtained by 
means of an ensemble average [13| accounting for the 
non- vanishing variance of the front position [lOj. Then 
the physical amplitude does no longer show geometric 
scaling [HEi. 
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APPENDIX A: MAPPING: FKPP^BK 

The fixed-coupling BK equation for N — N(Y, L) in 
the second-order ("diffusive") approximation: 



dN 

w 



d 2 L dN , 
X2 dI 2+Xl ^L +* oN ~ N ' ( A1 ) 



with: 



Xo = 7cV(7c)/2, 

Xi = 7cX"(7c) +X(7c)/7c 

Xo = x"(7c)/2, (A2) 

is mapped onto the FKPP equation for u = u(t,x): 

du d 2 u 9 , . „. 

m = o^ + u - u > (A3) 

by the relation: 

N{Y,L) = [ 7c 2 x"(7c)/2] x u[t(Y),x(Y,L)}, (A4) 

with: 

t = aY x [ 7c V( 7c )/2], 

x = aY x [ 7c 2 x"(7c) + X(7c)] + L x lc . (A5) 
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